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Abstract 



The intrinsically relativistic problem of a fermion subject to a pseudoscalar 
screened Coulomb plus a uniform background potential in two-dimensional 
space-time is mapped into a Sturm-Liouville. This mapping gives rise to an 
effective Morse-like potential and exact bounded solutions are found. It is 
shown that the uniform background potential determinates the number of 
bound-state solutions. The behaviour of the eigenenergies as well as of the 
upper and lower components of the Dirac spinor corresponding to bounded 
solutions is discussed in detail and some unusual results are revealed. An 
apparent paradox concerning the uncertainty principle is solved by recurring 
to the concepts of effective mass and effective Compton wavelength. 



1 Introduction 



The Coulomb potential of a point electric charge in a 1+1 dimension, con- 
sidered as the time component of a Lorentz vector, is linear (~ \x\) and so it 
provides a constant electric field always pointing to, or from, the point charge. 
This problem is related to the confinement of fermions in the Schwinger and 
in the massive Schwinger models 0-0 and in the Thir ring- Schwinger model 
0. It is frustrating that, due to the tunneling effect (Klein's paradox), there 
are no bound states for this kind of potential regardless of the strength of 
the potential 0-0. The linear potential, considered as a Lorentz scalar, is 
also related to the quarkonium model in one-plus-one dimensions 0-0. Re- 
cently it was incorrectly concluded that even in this case there is solely one 
bound state 0. Later, the proper solutions for this last problem were found 
0- [II]- However, it is well known from the quarkonium phenomenology in 
the real 3+1 dimensional world that the best fit for meson spectroscopy is 
found for a convenient mixture of vector and scalar potentials put by hand 
in the equations (see, e.g., |I2|). The same can be said about the treatment 
of the nuclear phenomena describing the influence of the nuclear medium 
on the nucleons ^3]-[22]- The mixed vector-scalar potential has also been 
analyzed in 1+1 dimensions for a linear potential as well as for a gen- 
eral potential which goes to infinity as \x\ — > oo |24| . In both of those last 
references it has been concluded that there is confinement if the scalar cou- 
pling is of sufficient intensity compared to the vector coupling. Although the 
vector Coulomb potential does not hold relativistic bound-state solutions, its 
screened version (~ e~^ x ^ x ) is a genuine binding potential and its solutions 
have been found for fermions |25j . The screened Coulomb potential has also 
been analyzed with a scalar coupling in the Dirac equation j2H] as well as in 
the Klein-Gordon equation with vector [27] and scalar [2E] couplings. 

The confinement of fermions by a pseudoscalar double-step potential [2H] 
and their scattering by a pseudoscalar step potential 01] have already been 
analyzed in the literature providing the opportunity to find some quite in- 
teresting results. Indeed, the two-dimensional version of the anomalous 
magnetic-like interaction linear in the radial coordinate, christened by Moshin- 
sky and Szczepaniak [3T] as Dirac oscillator, has also received attention. 
Nogami and Toyama [32], Toyama et al. [33J and Toyama and Nogami |34. 
studied the behaviour of wave packets under the influence of that conserving- 
parity potential whereas Szmytkowski and Gruchowski [33] proved the com- 
pleteness of the eigenfunctions. More recently Pacheco et al. [36J studied 
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some thermodynamics properties of the 1+1 dimensional Dirac oscillator, 
and a generalization of the Dirac oscillator for a negative coupling constant 
was presented in Ref. The two-dimensional generalized Dirac oscilla- 

tor plus an inversely linear potential has also been addressed jHH]- in recent 
papers, Villalba jHH] and McKeon and Van Leeuwen [40J considered a pseu- 
doscalar Coulomb potential (V = A/r) in 3+1 dimensions and concluded 
that there are no bounded solutions. The reason attributed in Ref. j3U| for 
the absence of bounded solutions it that the different parity eigenstates mix. 
Furthermore, the authors of Ref. jlHI assert that the absence of bound states 
in this system confuses the role of the it -meson in the binding of nucleons. 
Such an intriguing conclusion sets the stage for the analyses by other sorts 
of pseudoscalar potentials. A natural question to ask is if the absence of 
bounded solutions by a pseudoscalar Coulomb potential is a characteristic 
feature of the four- dimensional world. In Ref. [37] the Dirac equation in one- 
plus-one dimensions with the pseudoscalar power-law potential V = [i\x\ s 
was approached and there it was concluded that only for 5 > there can 
be a binding potential. Furthermore, we succeed in searching bound-state 
solutions for the pseudoscalar Coulomb potential (V = mu\x\) as well as for 
the pseudoscalar potential V = mux, even in the case uj < [3*7| . 

In the present paper we begin our presentation of the Dirac equation 
with the most general potential in 1+1 dimensions. This general approach 
allows us to show that the pseudoscalar potential has a clearly different be- 
haviour when compared with the scalar and vector potentials, even in the 
nonrelativistic limit of the theory. It was shown in Ref. J3S] that the prob- 
lem of a fermion under the influence of a pseudoscalar potential can always 
be mapped into a Sturm-Liouville problem for the upper and lower com- 
ponents of the Dirac spinor. Following that same methodology the present 
paper shows that the screened Coulomb potential plus a uniform background 
potential gives rise to an effective Morse-like potential, in the same way as 
the generalized Dirac oscillator plus an inversely linear potential gives rise 
to an effective quadratic plus inversely quadratic potential [38 . The exact 
bounded solutions are found, such as in Ref. j3H] , by transforming the Sturm- 
Liouville problem into Kummer's equation and expressing the upper and the 
lower components of the Dirac spinor in terms of confluent hypergeometric 
functions. In addition to their intrinsic importance as a new solution of the 
Dirac equation, this problem highlights the essential role of the uniform back- 
ground potential in furnishing bounded solutions and it might be relevant to 
studies of binding of neutral fermions by electric fields. Furthermore, the 
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result renders a new contrast to the result found in Ref . j3U] • I n the nonrel- 
ativistic quantum mechanics, the problem of a particle subject to the Morse 
potential (V — Vq [1 — exp(— ar)] 2 ) has been used to describe the vibrations 
of nuclei in homonuclear diatomic molecules |H]-jl21- This sort of problem 
for S-wave states is transformed into the problem of solving a transcendent 
equation which is only approximately solvable Nevertheless, the 

one-dimensional asymmetric Morse potential (—00 < r < 00) is an exactly 
solvable problem in the nonrelativistic quantum mechanics jlS]-jinj, even if 
its parameters are complex numbers |47j . 



2 The Dirac equation in a 1 + 1 dimension 

Thee 1+1 dimensional time- independent Dirac equation for a fermion of rest 
mass m reads 

= E^, H = cap + f3mc 2 + V (1) 

where E is the energy of the fermion, c is the velocity of light and p is the 
momentum operator. We use a = <7\ and (3 = 03, where o\ and 03 are Pauli 
matrices. For the potential matrix we consider 

V = W t + (3V S + aV e + Pj% (2) 

where 1 stands for the 2x2 identity matrix and fij 5 = 02- This is the most 
general combination of Lorentz structures for the potential matrix because 
there are only four linearly independent 2x2 matrices. The subscripts for the 
terms of potential denote their properties under a Lorentz transformation: t 
and e for the time and space components of the 2-vector potential, s and p 
for the scalar and pseudoscalar terms, respectively. It is worth to note that 
the Dirac equation is covariant under x — > —x if V e (x) and V p (x) change 
sign whereas V t (x) and V s (x) remain the same. This is because the parity 
operator P = exp(ir])P <j 3 , where 77 is a constant phase and Pq changes x 
into —x, changes sign of a and fi^ 5 but not of 1 and j3. 
Defining 

= exp Qa) K{x) = f dx ,] ^- (3) 
the space component of the vector potential is gauged away 

p + * = exp pip (4) 
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so that the time-independent Dirac equation can be rewritten as follows: 



Hip = EiP, H — a lC p + a 2 V p + a 3 (mc 2 + V s ) + W t (5) 

showing that the space component of a vector potential only contributes to 
change the spinors by a local phase factor. 

Provided that the spinor is written in terms of the upper and the lower 
components 

*=(t + ) («) 



the Dirac equation decomposes into: 



V t -E±V s ± mc 2 ) V>± = ihcxf)'^ ± iV p %l) T (7) 

where the prime denotes differentiation with respect to x. In terms of ip + 
and ip- the spinor is normalized as dx(\ip + \ 2 + \tp-\ 2 ) = 1 so that ip + 
and if)- are square integrable functions. It is clear from the pair of coupled 
first-order differential equations (J7J) that both ip + and have definite and 
opposite parities if the Dirac equation is covariant under x — ► —x. 

In the nonrelativistic approximation (potential energies small compared 
to mc 2 and E ~ mc 2 ) Eq. (J7J) becomes 

/ t2 t2 T/2 hV'\ 

Eq. (jHJ) shows that ?/>_ is of order v/c « 1 relative to tp + and Eq. © shows 
that ip + obeys the Schrodinger equation without distinguishing the contri- 
butions of vector and scalar potentials (at this point the author digress to 
make his apologies for mentioning in former papers Q29J- 30J) that the pseu- 
doscalar potential does not present any contributions in the nonrelativistic 
limit). 

It is remarkable that the Dirac equation with a nonvector potential, or a 
vector potential contaminated with some scalar or pseudoscalar coupling, is 
not invariant under V — ► F+const, this is so because only the vector potential 
couples to the positive-energies in the same way it couples to the negative- 
ones, whereas nonvector contaminants couple to the mass of the fermion. 
Therefore, if there is any nonvector coupling the absolute values of the energy 
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will have physical significance and the freedom to choose a zero-energy will 
be lost. This last statement remains truthfully in the nonrelativistic limit 
if one considers that such a contaminant is a pseudoscalar potential. It 
is also noticeable that the pseudoscalar coupling results in the Schrodinger 
equation with an effective potential in the nonrelativistic limit, and not with 
the original potential itself. Indeed, this is the side effect which in a 3+1 
dimensional space-time makes the linear potential to manifest itself as a 
harmonic oscillator plus a strong spin-orbit coupling in the nonrelativistic 
limit [3T]. The form in which the original potential appears in the effective 
potential, the V 2 term, allows us to infer that even a potential unbounded 
from below could be a confining potential. This phenomenon is inconceivable 
if one starts with the original potential in the nonrelativistic equation. 

As we will see explicitly in the next section, a constant added to the 
pseudoscalar Coulomb potential is undoubtedly physically relevant. As a 
matter of fact, it plays a crucial role to ensure the existence of bounded 
solutions. Nevertheless, the resulting potential does not present any nonrel- 
ativistic limit. Oddly enough, E « mc 2 only for the positive-ground-state 
solution in the strong-coupling regime of the theory. 



3 The pseudoscalar screened Coulomb poten- 
tial 

From now on we shall restrict our discussion to a pure pseudoscalar potential. 
For E 7^ imc 2 , the coupling between the upper and the lower components 
of the Dirac spinor can be formally eliminated when Eq. (JJj) is written as 
second-order differential equations: 



^ » - r[±] 



2m 

where 



— r± + v$ v± = E eff ^ (io) 



E 2 — m 2 c 4 

Eeff = 2mc 2 (U) 

V [±] - ^—±—V (12) 
K " ~ 2mc 2± 2mc 



These last results show that the solution for this class of problem consists 
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in searching for bounded solutions for two Schrodinger equations. It should 
not be forgotten, though, that the equations for tp + or tp_ are not indeed in- 
dependent because the effective eigenvalue, E e ff, appears in both equations. 
Therefore, one has to search for bound-state solutions for v}jj and V^jJ with 
a common eigenvalue. The Dirac eigenvalues are obtained by inserting the 
effective eigenvalues in (JTTJ). 

The solutions for E = ±mc 2 , excluded from the Sturm-Liouville problem, 
can be obtained directly from the Dirac equation (|7J). One can observe that 
such a sort of isolated solutions can be written as 



ip T = N T exp [=Ft>(x)] 

(13) 

2?7iC 

where N + and iV_ are normalization constants and v(x) = [ x dyV(y) /{he). 
One can check that it is impossible to have both components different from 
zero simultaneously on the same side of the x-axis. Of course normalizable 
eigenstates are possible only if v(x) has a distinctive leading asymptotic be- 
haviour. 

Now let us focus our attention on a pseudoscalar screened Coulomb po- 
tential in the form 

Tic 

g - g exp 



V 2A 




(14) 



where the coupling constants, go and g, are real numbers and A is a positive 
parameter related to the range of the interaction. Although g < gives rise 
to an ubiquitous repulsive potential in a nonrelativistic theory, the possibil- 
ity of such a sort of potential to bind fermions is already noticeable in the 
nonrelativistic limit of the Dirac theory (see Eq. Q). The presence of the 
uniform background potential, hcgo/(2X), is a distinguishing feature from 
the screened Coulomb potential in Refs. j2i]-j2H]- As one will see, it can not 
vanish for a binding pseudoscalar potential, and there is a minimum size for 
the coupling constant g before any bounded solution can be obtained, viz. 
| gfo | = 1- Evidently, the bound-state solutions for this potential present no 
nonrelativistic limit, hence we have to cope with an intrinsically relativistic 
bound-state problem. 



6 



We begin with the isolated solutions. In this case 



where e(x) stands for the sign function (e(x) = x/\x\ for x 7^ 0). It is not 
difficult to see that ip + and ^_ are not continuous at x — 0, and for that 
reason those solutions should be discarded. 

As for E 7^ 0, the effective potential becomes the Morse-like potential 

V$] = Aexp (-2^ + B ± exp + C (16) 

where 

2 2 2 2 2 

A = ^~>^ B ± = -£jL[g oT e(x)}, C = ^>0 (17) 
8X z m AX z m 8\ z m 

Before proceeding, it is useful to make some qualitative arguments re- 
garding the Morse-like potential and its possible solutions. The effective 
potential is able to bind fermions on the condition that B± < 0, contrariwise 
v}fj, or v}f}, or both of them would be repulsive everywhere. It follows 
that e(g)go > 1. It also follows that V^J < E e ff < C, where V^fj is the 
lowest value of v}fj. This implies that the Dirac eigenvalues corresponding 
to bounded solutions are in the range m 2 c 4 < E 2 < m 2 c 4 + 2mc 2 C. There 
is a spectral gap in the range E 2 < m 2 c 4 , and the energies belonging to 
E 2 > m 2 c 4 + 2mc 2 C correspond to the continuum. 

Note that the range of discrete Dirac eigenvalues enlarges as A decreases 
or | po| increases. As a matter of fact, when A decreases the potential well 
becomes deeper and beyond this it becomes narrower in such a way that 
A 2 !/^ is approximately constant, thus the capacity of the potential well to 
hold bound-state solutions remains almost the same. As for \go\, its increasing 
makes the potential well deeper almost without modifying its width. In this 
way, one expects that the parameter go determinates the number of allowed 
bounded solutions. 

The Morse-like potential presents, for B± < 0, a structure of asymmet- 
ric potential wells and the highest well governs the value of the zero-point 
energy. When g — > 0, v}fj has a single-well structure and nondegenerate 
Dirac eigenvalues are expected. When g — > oo, though, there appears a 



7 



double-well structure with an infinitely high and broad potential wall at the 
origin separating the two wells. In this last circumstance one would expect 
doubly-degenerate energy levels. 

Note that the parameters of the effective potential are related in such a 
manner that the change x — > —x induces the change V^fj — > v}jj without 
affecting the effective energy, signifying that \ip±{— x)\ behaves like |^ T (x)| 
and permitting us to focus our attention for a while on the positive side of 
the x-axis. 

Now let us do it quantitatively. Defining the dimensionless quantities z, 
/4ti) and ^ 



z = \g\ exp 



x 



hc\ \ 2A 



\ 2 



+ m 2 c 4 - E 2 



and using (fTT1|) - (fTT|) and (fT?)|l - fT7jl . one obtains the equation on the positive 
side of the x-axis 



^ + V4 + (-f-T + ^^ ± = (19) 



Now the prime denotes differentiation with respect to z. Following the steps 
of Refs. |2Zj and [2Hj, we make the transformation ip± = z~ l l 2 Q± to obtain 
the Whittaker equation [48J: 

4 ' z z 



$^ + ( -- + ^ + ' 2 | $± = (20) 



whose solution vanishing at the infinity becomes 

*± = N ± z^e^M (algj, b, z) (21) 

where N + and AL are constants, and M is the regular solution of the con- 
fluent hypergeometric equation (Kummer's equation) 
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zM" + {b-z)M' - 4(g) M = (22) 

with 

^l) = ^+\-^-(9oTl), b = 2u + l (23) 

Now we are ready to write the physically acceptable solutions on both sides 
of the x-axis by recurring to the symmetry x)\ ~ |^ T (x)| mentioned 

before. They are 



z v e~ z ' 2 



z»e- z ' 2 



9(-x)C^M (a", 6, z) + 9(+x)C^M (ag )s b, 
9(-x)D^M (ogj, 6, z) + 0(+z)L>WM (a^, 6, 



(24) 



where C^' and Z)M are constants, and 9(x) is the Heaviside function. 
The continuity of the wavefunctions (|24|) at x = furnishes 



(25) 



£>(+) £>I 



M ( a '|J)- g (^)^^o 



(26) 



where z = \g\. Substituting the solutions (}2"4")) into the Dirac equation (0) 
and making use of the recurrence formulas jlHj 



M'(a, b, z) = -Mia + 1, b + 1, z) 
b 



bM(a, 6, z) - bM(a - 1, b, z) - zM(a, b + l,z) = 

one has SIS db result 



(27) 



CM 



± 



A -E =F e(g)mc 2 



±e(g) 



(28) 



Together, (pfij) and (J2EJ lead to the quantization condition 
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> 1, and it represents a limit im- 



> m 2 c 4 , as has been anticipated by the preced- 



This last equation only holds for go/a[ 
posed on E 2 \ namely, E 2 
ing qualitative arguments. The transcendental eigenvalue equation ({2Hj) is 
invariant under the transformation E — ► —E as well as under the com- 
bined transformations g$ — > —go and g — ► — and it can be solved easily 
with a symbolic algebra program by searching for solutions in the range 
m 2 c 4 < E 2 < m 2 c 4 + (hcg /2X) 2 , with e(g)g > 1. Fig. [T] shows the be- 
haviour of the positive-eigenenergies as a function of g for A equal to the 
Compton wavelength, i.e. A = h/(mc), and go = 5. There is a finite number 
of eigenvalues. Notice that for very small values of g, at least one energy 
level emerging from the continuum comes into existence. The ground state 
tends asymptotically to ±mc 2 for arbitrary large g whereas the excited states 
tend to be closely bunched in pairs with eigenvalues independent of g. Due 
to the appearance of the potential barrier for large g, this sort of two- fold 
degeneracy comes as no surprise. 

Now we return our attention to (|24j) . where it remains to evaluate the 
constants CM and £>M. Using (J23J) and (J2HJ) one obtains 



Hg)] 



-e{g) 



ihc 

A E + e(g)mc 2 



he 



Mg)] 
l -e(g) 



A ^[E 2 



m 2 c A 



e{g) 



(30) 



h = m/P-^ [+] 

V -£/ + mcr 

and the constant C' + l is to be fixed by the normalization condition. It is 
noticeable in (j2H a somewhat left-right symmetry involving ip + and ■?/>_, 
such a symmetry is not exact inasmuch as C^' 7^ The upper and lower 

components of the spinor have the same number of zeros, or nodes, and as 
a consequence of such a quasi-symmetry the zeros of ip + and ^_ exhibit, 
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of course, an exact left-right symmetry. In whatever manner Eq. (|3U|) of 
course, is invariant under the change E — > — E provided \C^\ «-> |-D^|, 
hence one can conclude that |-0 + (±a;)| <-> \ip-(^fx)\ in such a way that the 
position probability density on the right side of the x-axis transforms into the 
position probability density on the left side of the x-axis, and vice versa. As 
for the combined transformations go — ► —go and g — ► —g, one has \C^\ <-> 
|£> [±] | with the proviso that m — ► — m. Recalling that the existence for 
bounded solutions demands that e{g)go > 1, this last set of mathematical 
transformations allow us to concentrate our attention on the case g > 0. Figs. 
I2MU illustrate the behaviour of the upper and lower components of the Dirac 
spinor, IV'+I 2 and \ip-\ 2 , and the position probability density, |^| 2 = l^+l 2 + 
IV'-I 2 , as a function of x, for the positive-energy solutions with g = 2, g = 5 
and A = h/(mc). Figs. do the same with g = 14. The normalization 
constant, CM, was fixed by numerical computation. Comparison of these 
figures shows that is larger than (for i? > 0) and that the fermion 
tends to avoid the origin as g increases. Furthermore, tends to concentrate 
at the left (right) region when E and g have equal (different) signs. This 
happens due to the discontinuity of the effective Morse-like potential at x — 0, 
given by v}fj(0+) — v}fj(0-) = ± h 2 g / (2X 2 m) . Note from all these figures 
that the quantum number n qualifies the number of nodes of ip+ an d ip-- 
Note also that the quasi-degeneracy of the second- and third-excited states 
for g = 14 reveals a near equality of the position probability densities for 
those states, as should be expected. 

A numerical calculation of the uncertainty in the position for the ground- 
state solution (with m — c — h — 1 and go = 5) furnishes 0.736 and 0.469 
for g = 2 and 14, respectively. Here we have purposely shown an odd fact 
for g = 14. It seems that the uncertainty principle dies away provided such 
a principle implies that it is impossible to localize a particle into a region 
of space less than half of its Compton wavelength (see, e.g., Ref. [49J). 
This apparent contradiction can be remedied by recurring to the concepts of 
effective mass and effective Compton wavelength. Indeed, the third line of 
()18|) suggests that we can define the effective mass as 



m 



eff 



r 2 +(f£) <-> 



Hence, the effective Compton wavelength can be defined as A e // = h/(m e ffc) 
so that the minimum uncertainty consonant with the uncertainty principle 
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is given by X e ff/2. Therefore, the uncertainty in the position can shrink 
without limit as \go\ increases or A decreases. It means that the localization 
of the fermion does not require any minimum value as | go \ — > oo or A — > in 
order to ensure the single-particle interpretation of the Dirac equation. 

4 Conclusions 

We have succeed in searching for Dirac bounded solutions for neutral fermions 
by considering a pseudoscalar screened Coulomb potential in 1+1 dimensions. 
The satisfactory completion of this task has been alleviated by the methodol- 
ogy of effective potentials which has transmuted the question for E ^ ±mc 2 
into Schrodinger-like equations with effective Morse-like potentials. There 
are no isolated solutions (E = ±mc 2 ) for this sort of potential. Nevertheless, 
as the coupling constant g becomes extremely large, i.e. \g\ ^> 1, the energy 
levels corresponding to the ground-state solution end up close to ±mc 2 . 

The uniform background potential, beyond the screened Coulomb poten- 
tial, is a sine qua non condition for the existence of bound-state solutions. 
Curiously enough, it plays an essential role not only to determinate the at- 
tractiveness of the effective Morse-like potential but also for establishing the 
band of allowed discrete Dirac eigenvalues. Furthermore, due to the fact 
that there is no atmosphere for the production of fermion-antifermion pairs, 
a neutron fermion embedded in this uniform background acquires an effective 
mass which permits that it can be strictly localized. 

In addition to their intrinsic importance as new solution of a fundamental 
equation in physics, the solutions obtained in this paper might be of relevance 
to the confinement of neutral fermions in a four- dimensional world. Further- 
more, they render a contrast to the result and conclusions found in j3U| : there 
are bound-state solutions for neutral fermions interacting by a pseudoscalar 
screened Coulomb potential in 1+1 dimensions, notwithstanding the spinor 
is not an eigenfunction of the parity operator. 
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Figure 2: \ip + \ 2 (full thin line), |t/>_| 2 (dashed line) and \ip\ 2 = \ip + \ 2 + \ip-\ 2 
(full thick line) as a function of x, corresponding to the positive-ground-state 
energy (n = 0) with X — h/ (mc), g — 2 and g — 5 (m = c = = 1). 
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Figure 3: The same as in Fig. 2 corresponding to positive- first-excited-state 
energy (n = 1). 
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Figure 4: The same as in Fig. 2 corresponding to positive-second-excited- 
state energy (n = 2). 



19 




20 




Figure 6: The same as in Fig. 2 for n = 1 and g = 14. 
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Figure 7: The same as in Fig. 2 for n = 2 and g = 14. 
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